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Abstract 

There is currently a great deal of interest in the invisibility cloaks recently proposed by 
Pendry et al. that are based in the transformation approach. They obtained their results 
using first order transformations. In recent papers Hendi et al. and Cai et al. considered 
invisibility cloaks with high order transformations. In this paper we study high order electro- 
magnetic invisibility cloaks in transformation media obtained by high order transformations 
from general anisotropic media. We consider the case where there is a finite number of spher- 
ical cloaks located in different points in space. We prove that for any incident plane wave, 
at any frequency, the scattered wave is identically zero. We also consider the scattering of 
finite energy wave packets. We prove that the scattering matrix is the identity, i.e., that 
for any incoming wave packet the outgoing wave packet is the same as the incoming one. 
This proves that the invisibility cloaks can not be detected in any scattering experiment 
with electromagnetic waves in high order transformation media, and in particular in the first 
order transformation media of Pendry et al. We also prove that the high order invisibility 
cloaks, as well as the first order ones, cloak passive and active devices. The cloaked objects 
completely decouple from the exterior. Actually, the cloaking outside is independent of what 
is inside the cloaked objects. The electromagnetic waves inside the cloaked objects can not 
leave the concealed regions and viceversa, the electromagnetic waves outside the cloaked 
objects can not go inside the concealed regions. 
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As we prove our results for media that are obtained by transformation from general 
anisotropic materials, we prove that it is possible to cloak objects inside general crystals. 

1 Introduction 

Recently fi\ considered electromagnetic invisibility cloaks based in the transformation method 
that offer the theoretical and practical possibility of hiding objects from observation by elec- 
tromagnetic waves. The results in [1] were obtained in the geometrical optics approximation. 
Numerical simulations were reported in [2] and [3] and a experimental verification of cloaking 
was given in P] . In the case of one spherically symmetric cloak [5] proved by an explicit cal- 
culation in spherical coordinates that the scattered wave is identically zero at each frequency 
and that the electromagnetic waves outside the spherical cloaked object can not enter the 
concealed region. These papers considered first order transformations. The recent papers 
[6l[7] consider high order transformations that offer new possibilities. In all these studies the 
transformation media were obtained by transformations from isotropic media. 

In this paper we consider electromagnetic invisibility cloaks in high order transformation 
media that are obtained by high order transformations from general anisotropic media. 
Moreover, we assume that there are several cloaked objects located in different points in 
space. 

We prove that for any incident plane wave, at any frequency , the scattered wave is 
identically zero. This generalizes, with a different proof, the results obtained in Note 
that in our case separation of variables can not be used as there is no symmetry since we 
have a finite number of cloaks in different points in space and also because we transform 
from general anisotropic media. 

We also consider the scattering of finite energy wave packets. We prove that the scattering 
matrix is the identity, i.e., that for any incoming wave packet the outgoing wave packet is 
the same as the incoming one. 

These results prove that the invisibility cloaks can not be detected in any scattering ex- 
periment with electromagnetic waves in high order transformation media, and in particular, 
in the first order transformation media of [1]. 
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We also prove that the high order invisibihty cloaks, as well as the first order ones, 
cloak passive and active devices. The cloaked objects completely decouple from the exte- 
rior. Actually, the cloaking outside is independent of what is inside the cloaked objects. 
The electromagnetic waves inside the cloaked objects can not leave the concealed regions 
and viceversa, the electromagnetic waves outside the cloaked objects can not go inside the 
concealed regions. 

In [S] we considered first order transformation media. Here we present the generalization 
to high order transformations, and we add new results. 

The fact that the electromagnetic waves inside the cloaked objects can not leave the 
concealed region, that we had already proved in [8] in the case of first order transformations, 
has recently been verified in the case of one spherical cloak in a first order transformation 
medium, by an explicit computation in spherical coordinates, by [9J, where also a physical 
mechanism based in surface voltages is presented to physically explain why the electromag- 
netic waves cannot leave the concealed region. 

Our results are based in Von Neumann's method of self-adjoint extensions. This is a 
very powerful technique that allows us to settle in an unambiguous way the mathematical 
problems posed by the singularities of the inverse of the permittivity and the permeability 
of the transformation media in the boundary of the cloaked objects. It also allows us to 
identify the appropriate boundary condition when cloaking is formulated as a boundary value 
problem. Namely, that the tangential components of the electric and the magnetic fields have 
to vanish at the outside of the boundary of the cloaked objects. See Remark 12.41 We have 
already proven this result for first order transformations in |8|. This boundary condition is 
self-adjoint in our case because the permittivity and the permeability are degenerate at the 
boundary of the cloaked objects. 

As we prove our results for media that are obtained by transformation from general 
anisotropic materials, we prove that it is possible to cloak objects inside general crystals. 

As it is often the case in the papers on electromagnetic invisibility cloaks, I make the 
assumption that the media are not dispersive. This is a widely used idealization. As is well 
known, metamaterials are dispersive, and, furthermore, when the the permittivity and the 
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permeability have eigenvalues less than one, dispersion comes into play in order that the 
group velocity does not exceeds the speed of light. This idealization means that we have to 
take a narrow enough range of frequencies in order that we can analyse the cloaking effect 
without taking dispersion into account. In practice this means that cloaking will only be 
approximate. 

For a related method for cloaking in two dimensions see [lOj. For earlier results in cloaking 
for conductivity problems see [TT] , the references quoted there, and fTZl [T3] . In [T3j cloaking 
is studied in the context of the Dirichlet to Neumann map. For a cylindrical invisibility 
cloak with first order transformation see {T5\. See also UM and ITTI for other related results 

The paper is organized as follows. In Section 2 we introduce our formalism, we give 
the definition of solutions with locally finite energy and we obtain the cloaking boundary 
condition. Moreover, we prove that the electromagnetic waves inside the cloaked objects can 
not go outside and viceversa. In Section 3 we prove that the scattered waves are zero for all 
frequencies and all incoming plane waves. In Section 4 we prove that the scattering matrix 
is the identity for all incoming wave packets. In Section 5 we give the proof of Theorems 
12. 3[ and 12.51 In Section 6 we discuss generalizations of our results. Finally, we give a brief 
conclusion and outlook where we also comment on cloaking objects inside general anisotropic 
media, in particular inside general crystals. 

2 Electromagnetic Cloaking 

Let us consider Maxwell's equations in M^, in the time domain. 



V X E 



-— B, V X H = — D, 

dt ' dt ' 



(2.1) 



V B 



0,V-D 







(2.2) 
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and in the frequency domain, assuming a periodic time dependence of E, H given by 
e~™*, with uj the frequency, 

V X E = iujB, V X H = -iuj'D, 00 ^ 0, (2.3) 
V-B = 0,V-D = 0, (2.4) 

where we have suppressed the factor e~*^* in both sides. Note that (12.41) follows from (12.31) . 

We study the propagation of electromagnetic waves -that satisfy Maxwell's equation- in 
the case where there is a finite number of high-order spherical invisibility cloaks located in 
different points in space. 

For simplicity let us first consider the case where there is only one cloak located at 
X = 0. See Figure [TJ We designate the Cartesian coordinates of x by x'^, A = 1, 2, 3. To 
define the high-order transformation media we introduce another copy of M.^, denoted by 
Rq. The points in Rq are denoted by y with coordinates y^,X = 1,2,3. We designate by 
X = x/|x|,y = y/|y| unit vectors. Consider the following transformation from Rq \ {0} to 

r3 mil], 

x = x(y) = /(y):=^7(|y|)y. (2.5) 

In spherical coordinates this transformation changes the radial coordinate but leaves the 
angular coordinates constant, i.e., |x| = (7(|y|),x = y. Given < a < 6 we wish that this 
transformation sends the punctuated ball < |y| < b onto the concentric shell a < |x| < 6, 
that it is the identity for |y| > 6 and that it is one-to-one. Then, we assume that g satisfies 
the following conditions. 

DEFINITION 2.1. For any positive numbers a, b with < a < b, we say the g is a cloaking 
function in [0,6] if g{p) is twice continuously differentiable on [0,6], (7(0) = a,g{b) = b, and 
g'ip) := A^(p)>0,pG[0,6]. 

These high order transformations were introduced in [7] for the case of a cylindrical cloak. 
They imposed the condition g'{b) = 1. Since for our spherical cloaks we do not need this 
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condition we do not assume it. We define, 



X = x(y) = /(y) := ^(|y|)y, forO < |y| < b, 
X = x(y) := y, for |y| > b. 



(2.6) 



Witli tliese conditions (12.61) is a bijection from Rq \ {0} onto \ -Ba(O), where by 



we denote the closed ball of center xq and radius r. Moreover, it blows up the point onto 
the sphere |x| = a. It sends the punctuated ball < |y| < 6 onto the concentric shell 
a < < b and it is the identity for |y| > b. It is twice continuously differentiable away 
from the sphere |y| = b, where it can have discontinuities in the derivatives depending on 
the values of the derivatives of g at b. 

In [7] the quadratic case 



with p G M was discussed in connection with a cylindrical cloak in an approximate transfor- 



The closed ball := {x G : |x| < a} is the region that we wish to conceal, and we 
call it the cloaked object. The spherical shell a < |x| < 6 is the cloaking layer. The union of 
the cloaked object and the cloaking layer is the spherical cloak, that in this case is just the 
closed ball of center zero and radius b. The domain |x| > 6 is the exterior of the spherical 
cloak. 

To have several cloaks we just put a finite number of these spherical cloaks in differ- 
ent points in space at a finite positive distance from each other, in order that they do 
not intersect. See Figure [H Let us take as centers of the cloaks points Cj G IR^,J = 
1,2,- --A^ where is the number of cloaks and Cj ^ C;,j 1,1 < j,l < N. We take 
< aj < bj, and cloaking functions (yfj that satisfy the conditions of Definition 12.11 for 
ttj, bj,j = 1,2,3 ■■■ N, and we define the following transformation from Mq \ {ci, C2, ■ ■ ■ , cjy} 



Br{xo) := {x G : I a; - Xo| < r 



} 



(2.7) 



dip) = T+Pip-b) P + a 



mation medium. In [I] the first order case g{p) = ^-ff^p + a was considered. 
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to R^. 



X = x(y) = /(y) := Cj + gj{\y - Cj\) y - Cj, y G fi;,,(Cj), j = 1, 2, ■ ■ ■ , A^, 

(2.8) 

X = x(y) = /(y) := y,y e M-j \ Uf^,B,^{c^), 
where Bb^{cj) are balls in Rq. 

The cloaked objects that we wish to conceal are given by, 

Kj := {xG : |x-Cj| < aj} ,j = 1,2,--- ,N. (2.9) 

The concentric spherical shells aj < |x — Cj| < bj,j = 1,2, ■■ ■ , are the cloaking layers. 
The spherical cloaks are the balls Bi,-[cj) in . We denote by K the union of all the cloaked 
objects, 

K := uf^^K,. (2.10) 

The domain 

Mnuf=ii?.,(c,) (2.11) 

is the exterior of the all the spherical cloaks. We assume that the spherical cloaks are at a 
positive distance of each other, 

min distance (^Bb.{cj), Bb^ci)) > 0, j 7^ /, j, / = 1, 2, ■ ■ ■ , N. 

Denote 

Qo ■■= \ {ci, C2, ■ ■ ■ , c;v}, ft:=M.^\ K. 

Then, f l2.8p is a bijection from VLq onto fi, and for j = 1, 2, ■ ■ ■ , it blows up the point 
Cj onto the sphere |x — Cj| = Oj. It sends the punctuated ball < |y — Cj| < 6^ onto the 
concentric shell < |x — Cj| < hj and it is the identity for y G Mg \ interior {u^^^Bb-^Cj)^. 
It is twice continuously differentiable away from the spheres |y — Cj| = 6j, where it can have 
discontinuities in the derivatives depending on the values of the derivatives of gj at hj. 

For any open set O and for any n = 1, 2, ■ ■ ■ , let us denote by C"'(0) the set of all 
C— valued functions that are continuous together with all its derivatives of order up to n, 
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and by Cq{0) the functions in C"(0) that have compact support in O, i.e. such that the 
closure of the set of points where they are different from zero is bounded in O. In other words, 
the closure of set of points where the functions are different from zero is bounded, and they 
are zero in a neighborhood of the boundary of O. By C{0), Co(0), we denote, respectively, 
the continuous functions in O and the continuous functions with compact support in O. 

We denote the elements of the Jacobian matrix by Ay, 

Note that A^, e [Qq \ uf^^dBb^{cj)), and that it can have discontinuities on U^—idBf,^ (cj) 
depending on the derivatives of gj at bj. We designate by A^' the elements of the Jacobian 
of the inverse bijection, y = y(x) = /~^(x), 

eC^{Q\ Uf=i95,^,(cj)), and it can have discontinuities on UjLidBi,^{cj) depending on 
the derivatives of gj at bj. Let us denote by A the determinant of the Jacobian matrix (I2.12p . 
See (15. 2 p for the calculation of A in closed form. Note that it diverges at dK = dQ. 

We take here the material interpretation and we consider our transformation as a bijection 
between two different spaces, Qq and Q. However, our transformation can be considered, as 
well, as a change of coordinates in Qq. Of course, these two points of view are mathematically 
equivalent. This means, in particular, that under our transformation Maxwell's equations 
in Qo and in Q have the same invariance that they have under change of coordinates in 
three-space. See, for example, [18]. Let us denote by Eq, Hq, Bq, Dq, ^q*^, /Xo'^, respectively, 
the electric and magnetic fields, the magnetic induction, the electric displacement, and the 
permittivity and permeability of Qq. The, Sq'^^^q^, are positive, hermitian matrices that are 
constant in VLq. 

The electric field is a covariant vector that transforms as, 

EA(x) = A^'(y)Eo,A'(y). (2.14) 
The magnetic field H is a covariant pseudo- vector, but as we only consider space transfor- 
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mations with positive determinant, it also transforms as in fl2.14p . The magnetic induction 
B and the electric displacement D are contravariant vector densities of weight one that 
transform as 

B\^) = {A{y))-'Al,{y)B^'{y), (2.15) 

with the same transformation for D. The permittivity and permeability are contravariant 
tensor densities of weight one that transform as, 

e'^i^) = (A(y))-^ A^,(y) A::,(y)4'^'(y), (2.16) 

with the same transformation for /i'^'^. Maxwell's equations (12. 1112. 41) are the same in both 
spaces Q and Qq- Let us denote by e\u, fixu,soxu, fJ'OXu, respectively, the inverses of the 
corresponding permittivity and permeability. They are covariant tensor densities of weight 
minus one that transform as, 

e,,(x) = A(y)Af (y)<(y)5oAv(y), /iA.(x) = A(y)Af (y) ^(y) /.oAV'(y)- (2-17) 
Note that 

det e^'' = A-^ det det /i^^ = A^^ det jj^", (2.18) 

det exu = A det Eqxu, det nxu = A det Hoxu- (2.19) 

Then, by (12.181 12.19[ 15. 2p the matrices e^'^,fj.^'' are degenerate at dK and the matrices 
£xu,fJ'Xu are singular at dK. 

We face now the problem that as the e'^'' and /x'^'^ are degenerate at the boundary of the 
cloaked objects K we have to make precise what do we mean by a solution to Maxwell's 
equations. As we will see the standard rules that we apply in non-degenerate situations 
do not apply here. These type of problems are not unusual in mathematical physics and 
there is a well established method to deal with them. Namely, Von Neumann's theory of 
self-adjoint extensions fl9[ [20]. Let us first consider the problem in Q. We write Maxwell's 
equations in Schrodinger form. For this purpose we denote by e and /x, respectively, the 
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matrices with entries exu and fix^. Recall that (V x E)'^ = s^'^^-^Ep, where s^^'' is the 
permutation contravariant pseudo-density of weight —1 (see section 6 of chapter II of [IB], 
where a different notation is used). 

We define the following formal differential operator, 

aJl\^.{ ) . (2.20) 



H y \^ -;uV X E 
Here, as usual, we denote, eV x H := exu{^ x H)'', and yuV x E = fixuC^ x E)"". 
Equation (12. ip is equivalent to, 

4(^)=an(^). (2.21) 



dt\ii J V H 

and equation 02.31) is equivalent to 



^( Si =«nf S ). (2.22) 



H ; V H 

Note that since the matrices e, /i are singular at dQ the operator an has coefficients that are 
singular at dQ. This is the reason why we have to be careful when defining the solutions. 

We have to define equation (12.201) in an appropriate linear subspace of the Hilbert space 
of all finite energy fields in Q, that we define now. We designate by Hqe the Hilbert space 
of all measurable, C^— valued functions defined on Q that are square integrable with the 
weight e^" and the scalar product, 

(E«,E(2))^^:= / Ef^e^^^^Frfxl (2.23) 
Jq 

Moreover, we denote by Tinn the Hilbert space of all measurable, C^— valued functions 
defined on Q that are square integrable with the weight /i'*"' and the scalar product, 

(H«,H(2))^^ := / i/fV'^^rfxl (2.24) 
Jn 

The Hilbert space of finite energy fields in Q is the direct sum 

nn:=nnE®nnH. (2.25) 
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For any open set O the spaces C(0), Co(0), C"(0), Cq(0), = 1,2 - ■■ are defined as 
the spaces C(0), Co(0), C"(0), Co"(0), n = 1, 2 ■ ■ ■ but with C^- valued functions. 

We first define in a nice set of functions where it makes sense, that we take as C\{VL). 
In physical terms this means that we start with the minimal assumption that Maxwell's 
equations are satisfied in classical sense away from the boundary of Vt. with domain 
-D(an) := CQ(fi) is a symmetric operator in Hn, i.e. C a^. To construct a unitary 
dynamics that preserves energy we have to analyze the self-adjoint extensions of a^, what in 
physical terms means that we have to make precise in what sense Maxwell's equations are 
solved up to dfl. In other words, to construct finite-energy solutions of fl2.21l) . 



E 
H 



with constant energy 

(E(t), E(t))^^ + (H(t), H(t))^^ = (E(0), E(0))^^ + (H(0), H(0))^^ < oo, 

we have to demand that the initial finite energy fields , (E(0), H(O))"^ belong to the domain 
of one of the self-adjoint extensions of a^. The key issue is that has only one self-adjoint 
extension, i.e. it is essentially self-adjoint. Before we state this result in a precise way 
in Theorem 12.31 let us discuss its physical consequences. Let us denote by the unique 
self-adjoint extension. 

We denote by kernel the null subspace of Aq, i.e.. 



kerneMo := <! ( ^ ) G D{An) : ^ 



and by 



'Hn± ■= (kerneMo)" 



the orthogonal complement in Tin of kernel Aq. Equations (12.21) are satisfied for all times 
if and only if (E(0), H(O))'^ G (kerneMn) . Moreover, the unique finite energy solutions to 
(12. H 12.21) with constant energy are constructed as follows. 

We take any 

G nn± n D{An) (2.26) 



E 
H 
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and we obtain the finite energy solution to Maxwell's equations fl2.H 12. 2p as 




(2.27) 



This is the unique finite energy solution with constant energy, and with initial value at t = 
given by fl2.26p . Note that as e~**^"?in± C 7in± equations (12. 2p are satisfied for all times if 
they are satisfied at t = 0. The unitary group e~**"^" is defined via functional calculus, but 
we can think of it as just the operator that gives us the unique solution. We can consider 
more general solutions by means of the scale of spaces associated with Aq, but we do not go 
into this direction here. 



Solutions to (12. 3112. 4p in general do not have finite energy because they do not have 

enough decay at infinity to be square integrable over all Q. Then, we only require that 
they are of locally finite energy in the sense that the electric and the magnetic fields are 
square integrable over every bounded subset of Q, respectively, with the weight e'^'^, and 
/i^^. Moreover, in order that the problem (12. 3112. 4p is well-posed -in the sense that it is 
self-adjoint- the solutions with locally finite energy have to be locally in the domain of the 
only self-adjoint extension of a^, that is to say, they have to be in the domain of Aq when 
multiplied by any continuously differentiable function with support in a bounded subset of 
Q. Hence, we define, 

DEFINITION 2.2. (Solutions with Locally Finite Energy) We say that the fields (E, H)'^ 
are a solution to Ii2.3^2.4\ ) with locally finite energy in Q if they solve Ii2.3\2.4\ ) in distribution 
sense in Q, if, furthermore, for every bounded set O G Q 



and if for any continuously differentiable function with bounded support in Q, (E, H)-^ G 



Similarly, given any bounded set O G Q we say that the fields (E, H)"^ are a finite 
energy solution in O if (12.280 is satisfied, if the (E, H)^ are a solution to (12.31 12. 4p in 
distribution sense in O and if for any continuously differentiable function (p with support in 





(2.28) 



D{An). 



O, (P{E,Hy eD{An). 
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To state Theorem 12.31 and to further analyze our problem we consider now Maxwell's 
equations in Rq and we define the Hilbert spaces of electric and magnetic fields with finite 
energy. The Eq, Hq, Bq, Dq, were defined in ^o, but since Mq \ = {cj}^^^ is of measure 
zero, we can consider them as defined in Mq, what we do below. 

We denote by Hqe the Hilbert space of all measurable, square integrable, C^— valued 
functions defined on Mq with the scalar product, 

E«eo'^EjJdy3. (2.29) 

■0 

We similarly define the Hilbert space,?-^oH, of all measurable, square integrable, C'^— valued 
functions defined on Mq with the scalar product, 

(H«,Hf) = / ffSV-^^Tdy'. (2.30) 

^ ''OH J^3 

The Hilbert space of finite energy fields in Mq is the direct sum 




T^o '■= T^oE ® T^oH ■ (2.31) 

We now write Maxwell's equations in Mq in Schrodinger form. As before we denote by Eq 
and fiQ, respectively, the matrices with entries Eqxi, and fj^oxu- By oq we denote the following 
formal differential operator. 



Eo \ ■ f SqV X Ho 



Then, equation (12.11) in Mq is equivalent to. 



^ / Eo \ / Eo 



Let us denote by C[^(Mq) the set of all C^— valued continuously differentiable functions 
on Rq that have compact support in Mq . Then, ao with domain C[f (Mq) is a symmetric 
operator in Ho, i.e., oq C Oq. Moreover, it is essentially self-adjoint in Hq, i.e., it has only 
one self-adjoint extension, that we denote by Aq. Its domain is given by. 



^(^o) = <! ( Hq ) ^ Hq ) ^ ^4 ' (2.34) 
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and, 

A. ( Ho ) = ^0 ( Ho ) ' ( Ho ) ^ (2.35) 
where the derivatives in the right-hand sides of fl2.34[ I2.35P are taken in distribution sense in 
Mq. These results follow easily from the fact that -via the Fourier transform- is unitarily 
equivalent to multiplication by a matrix valued function that is symmetric with respect 
to the scalar product of Tio- Moreover, it follows from explicit computation that the only 
eigenvalue of Aq is zero, that it has infinite multiplicity, and that. 



no± := (kernel Aq)" = | J ^ T^o : ^o. = 0, ^/"o ^o. = 0| . (2.36) 

Furthermore, Aq has no singular-continuous spectrum and its absolutely-continuous spec- 
trum is M. See, for example, [2T], [22] . 

Taking any 

( Ho ) e^o±nD(Ao) (2.37) 
we obtain a finite energy solution to Maxwell's equations (12. ![ I2.2p in Mj] as follows 

Ho ) = ( Ho ) • ^'-''^ 

This is the unique finite energy solution with initial value at t = given by (12.371) . Note 
that as e~**^''7io± C 7io± equations (12. 2p are satisfied for all times if they are satisfied at 
t = 0. 

We denote by Ue the following unitary operator from Hqe onto Hue, 

(f/ijEo);, (x) := Afi?oA'(y), (2.39) 
and by Uh the unitary operator from Hqh onto Hqh, 



^A^OA'(y). 



(2.40) 



Then, 



U := Uf 



U 



H 



(2.41) 
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is a unitary operator from Hq onto Hq. 

We prove the following theorem in Section 5. 

THEOREM 2.3. The operator aQ is essentially self-adjoint, and its unique self-adjoint 
extension, Aq, satisfies 

An = UAoU*. (2.42) 

Furthermore, Aq has no singular- continuous spectrum and its absolutely- continuous spectrum 
is M. The only eigenvalue of Aq is zero and it has infinite multiplicity. Moreover, 



E^ ^-Av7Ti_n ^ ..Xu . 



nn± := (kernel Anr = | r J ^ " d^,' " ^ °' d^/ " °J ' ^^'^^^ 

The facts that is essentially self-adjoint and that its unique self-adjoint extension Aq is 
unitarily equivalent to the generator Aq of the homogeneous medium are strong statements. 
They mean that the only possible unitary dynamics in Q that preserves energy is given by 
(12.271) and that this dynamics is unitarily equivalent to the free dynamics in Mq given by 
(12.381) . In fact, dQ acts like a horizon for electromagnetic waves propagating in Q in the 
sense that the dynamics is uniquely defined without any need to consider the cloaked objects 
K = U^^^Kj. As we will prove below this implies, in particular, electromagnetic cloaking 
for all frequencies in the strong sense that the scattered wave is identically zero at each 
frequency and that the scattering operator in the time domain is the identity. 

Formulating cloaking as a boundary value problem in Q is of independent interest. For 
this purpose we introduce the following boundary condition. 

REMARK 2.4. ( The Cloaking Boundary Condition) 

Let (E, H)^ be a solution with locally finite energy in Q. According to Definition ^ 
they have to be locally in the domain of Aq. By (I2.42p this implies that. 



E ^ U(^^ ) (2.44) 



H y V Ho 

with (Eo,Ho)^ locally in the domain of Aq, i.e., (Eo,Ho)^ is in the domain of Aq when 
multiplied by any function in Co(Mq). It follows from (I2.44p . and (15. ip that the solutions 
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with locally finite energy have to satisfy the cloaking boundary condition, 

E X n = 0, H X n = 0, in dQ = dK+, (2.45) 

where dK^ is the outside of the boundary of the cloaked objects and n is the normal 
vector to dK^. 

Note that as Aq is the only self-adjoint extension of a^, this is the only possible self- 
adjoint boundary condition on 9/^+. It is self-adjoint because the matrices e,fi are singular 
at dK^. Then, cloaking as boundary value problem consists of finding a solution to (12. 3[ 
12. 4p in Q with locally-finite energy that satisfies the cloaking boundary condition given in 

dMSD. 

□ 

Let us now consider the propagation of electromagnetic waves in the cloaked objects. 
For this purpose we assume that in each Kj the permittivity and the permeability are 
given by e^" , , with inverses ejXu,^j\u and where ej,^j are the matrices with entries 
£j\u,fJ'jXu- Furthermore, we assume that < e^'^ (x.) , fi^'^ (x.) < C, x G Kj and that for 
any compact set Q contained in the interior of Kj there is a positive constant Cq such 
that det£:^'^(x) > Cq , det yU'^'' (x) > Cq,x G Q. In other words, we only allow for possible 
singularities of ej,fij on the boundary of Kj. 

We designate by TijE the Hilbert space of all measurable, C'^— valued functions defined 
on Kj that are square integrable with the weight Ej'^ and the scalar product, 

(^?'^f%,^-i/^-rW^- (2.46) 

Similarly, we denote by Ti-jH the Hilbert space of all measurable, C^— valued functions 
defined on Kj that are square integrable with the weight fij'' and the scalar product, 

(H?'.Hf),.„- X^^fj^r^^x'. (2.47) 



16 



The Hilbert space of finite energy fields in Kj is the direct sum 

Hj := HjE © njH, (2.48) 
and the Hilbert space in the cloaked objects K is the direct sum, 



The complete Hilbert space of finite energy fields including the cloaked objects is, 



H:=nn®HK. (2.49) 

We now write (12. ip as a Schrodinger equation in each Kj as before. We define the 
following formal differential operator, 

•^A^^-^ l^^'^^- (2.50) 



Equation (12. ip in Kj is equivalent to 

9^E, \_„ /E 



o o 

Let us denote the interior of Kj by Kj-= Kj \ dKj. Then, aj with domain Cl{Kj) is a 
symmetric operator in Tij. We denote. 



a := an © ax, where a^ := ©j=iaj, (2.52) 

with domain. 

The operator a is symmetric in ?-^. The possible unitary dynamics that preserve energy for 
the whole system, including the cloaked objects, K, are given by the self-adjoint extensions 
of a. We have that. 
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THEOREM 2.5. Every self-adjoint extension, A, of a is the direct sum of Aq and of some 
self-adjoint extension, A^, of a^, i-e., 

A = An®AK. (2.54) 

This theorem tells us that the cloaked objects K and the exterior Q are completely 
decoupled and that we are free to choose any boundary condition inside the cloaked objects K 
that makes ax self-adjoint without disturbing the cloaking effect in Q. Boundary conditions 
that make Ak self-adjoint are well known. See for example, [23], [21], [25] and [26]. 

It follows from explicit computation that zero is an eigenvalue of every Ax with infinite 
multiplicity and that. 



nK± := (kerneM;^)" ^ ^ j ^ : q^^^rEu = 0, Q^^f^RHu = j , (2.55) 

where we denote by e^(x) := ^^'^(x) for x G Kj, and /i;^(x) := /ij'^(x) for x G Kj,j = 
1, 2, ■ ■ ■ ,N. It follows that zero is an eigenvalue of A with infinite multiplicity and that, 

:= (kernel A)^ = Hn^ © (2.56) 

For any ip = ipn ® e n± D D{A), 



e-^*^<^ = e-'*^" ifin © e"'*^^ (2.57) 

is the unique solution of Maxwell's equations (12.11 12.21) with finite energy that is equal to (p 
at t = 0. This shows once again that the dynamics in Q and in K are completely decoupled. 
If at t = the electromagnetic fields are zero in Q, they remain equal to zero for all times, 
and viceversa. Actually, electromagnetic waves inside the cloaked objects are not allowed 
to leave them, and viceversa, electromagnetic waves outside can not go inside. In general, 
the solutions will be discontinuous at dK. This implies, in particular, that the presence of 
active devices inside the cloaked objects has no effect on the cloaking outside. In terms of 
boundary conditions this means that transmission conditions that link the electromagnetic 
fields inside and outside the cloaked objects are not allowed. Furthermore, choosing a par- 
ticular self-adjoint extension of the electromagnetic propagator of the cloaked objects can 
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be understood as choosing some boundary condition on the inside of the boundary of the 
cloaked objects. In other words, any possible unitary dynamics that preserves energy im- 
plies the existence of some self-adjoint extension in K, that can be understood as a boundary 
condition on the inside of the boundary of the cloaked objects. The particular self-adjoint 
extension, or boundary condition, that nature will take depends on the specific properties 
of the metamaterials used to build the transformation media as well as on the properties 
of the media inside the cloaked objects. Note that this does not mean that we have to put 
any physical surface, a lining, on the surface of the cloaked objects to enforce any particular 
boundary condition on the inside, since as wc already mentioned this plays no role in the 
cloaking outside. It would be, however, of theoretical interest to see what the self-adjoint 
extension in K, or the interior boundary condition, turns out to be for specific cloaked ob- 
jects and metamaterials. These results apply to the exact transformation media that we 
consider on this paper. However, the fact that there is a large class of self-adjoint extensions 
-or boundary conditions- that can be taken inside the cloaked objects could be useful in or- 
der to enhance cloaking in practice, where one has to consider approximate transformation 
media, as well as in the analysis of the stability of cloaking. 

3 Scattering at Fixed Frequency 

We now consider the scattering of plane waves by the cloaks. We have the following result. 

THEOREM 3.1. For any incident plane wave at any frequency the scattered wave is iden- 
tically zero. 

Proof: For simplicity, let us first consider the case where the 6^'^,!!^''' are isotropic, i.e., 
Eq'^ — e^5^'^,iJLQ^ — 11^5^'^ and that we have an incident field that propagates along the 
vertical axis, X3, with the electric field polarized along the Xi direction. 




where k — uj^s^ijP, and ei, 62 are unit vectors, respectively, in the and x"^ directions. 
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We have to find the solution with locally finite energy to (12.31 12.4p . (E, H)^, that outside 
the cloaks is equal to the sum of the incident plane wave and the scattered wave, i.e., 



E 
H 



X 



rpin \ / rpSC 

Hin )W+( H^c )(x),xG U^li5.,(c,). 



(3.1) 



Furthermore, the scattered wave (E^*", H^'")-^ has to satisfy the outgoing Silver- Miiller radi- 
ation condition. See, for example, [27]. We compute this solution in a simple way using the 
unitary equivalence of Aq and Aq. By (I2.42p the fields with locally finite energy. 



E 
H 



U 



(x), 



are a solution to (12.31 12.41) . Furthermore, since for x G \ U^^-Bj,^ (cj), A^, = 6y, we have 
that. 



E 
H 



LBL Q pi(kx^—ut) 

\ Y MO ^ 



xGMnuf=,i?,^.(c,), 



what proves that the scattered wave is identically zero. 

Let us now consider the case where e^^, are general anisotropic media. For a discussion 
of plane wave solutions in this case see, for example, [21], [22l EH]- A general incident plane 
wave is given by. 



jjin 



E^'^(k) e^i'^-^-'^'t) 
H^^(k) e*^''^"'^*) 



where k G M^, k 7^ and. 



k X E"^(k) ) = ufi^'H'J'ik), ( k X H"^(k) ) = -a;£^X''(k). 



We compute again the solution with locally finite energy using the unitary equivalence 
between An and Aq (12.421) . We have that. 



E \ , , / Ei'^(k)e^('^-y-"*) 

H H-(k)e^(— ) 
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are a solution to (12.31 12.41) with locally finite energy. As in the isotropic case, 



what proves that the scattered wave is also zero in the general anisotropic case. 

□ 

This theorem proves that we can not detect the cloaked objects K in any scattering 
experiment. 

4 Scattering of Wave Packets 

In this section we consider the scattering of finite energy wave packets. Of course, in practice 
one always sends a finite energy wave packet on the target, that has to have a narrow 
enough range of frequencies in order that we can neglect dispersion, as we mentioned in the 
Introduction. 

Let xn be the characteristic function of fi, i.e., Xn(x) = 1, x G fi, Xni?^) = 0, x G \ fi. 
We define, 



Eo \ \ / N . X / Eo 



H Ho j j ^"f"' ( Ho j 

By (15. 3p J is a bounded operator from Tio into Ti^. 
The wave operators are defined as follows, 

W± = s- lim e**^" Je-^*^°Po±, (4.2) 

provided that the strong limits exist, and where Po± denotes the projector onto 'Ho±- 

Let us designate by W^'^(]Rg) the Sobolev space of valued functions. We denote by / 
the identity operator on Hq. Then, 

LEMMA 4.1. The wave operators ^4-2^ exist and, 

W± = UPo±. (4.3) 
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Proof: Denote, 



W{t):= e**^" Je-^*^oPo±- 

By ([232D, for any e Ho 

W{t)ip = tPit) + UPo±^, (4.4) 

with 

^{t) := f/e**^° {U*J - I) e-^*^oPo±<^- (4.5) 

Since for |y| > R, with R large enough, our transformation, x = /(y), is the identity, 
X = y, and in consequence, ^^(y) = 6y for |y| > R, we have that, 

(U*J-I) = iU*J-I)xBnio)- (4.6) 

It follows that. 



s- lim ij{t) = Us- lim e**^°i?(t) (4.7) 

t— ►±00 t— >±oo 

with, 

m ■■= iU\J - I) XBnio)e-''^"Po±¥>. (4.8) 

We have that. 



(4.9) 



||^(t)||„„ < II JxB^(o)e-^*^°Po±</^||„ + \\xBnio)e-''''"Po±v\\no 

Then, as {Aq + i)-^Po± is bounded from Ho into W1'2(R3) [2l] [22], it follows from the 
Rellich local compactness theorem that 

XBnio) (Ao + iy^Poi. 
is a compact operator in Ho- Suppose that (p G D{Aq) Pi 7Yo±. Then, 

s- lim XB«(o)e-'*^°Po±<^ = s- lim Xb«(o)(^o + 0~'^o±e-'*^°(Ao + z)^ = 0, (4.10) 
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and whence, by (14. 9p . 

and it follows that in this case, 



s- lim i^{t) = 0, 

4— >itoo 



(4.11) 



s- lim ip{t) = 0. (4.12) 

f— >itoo 

By continuity this is also true for ip G Ho±. 
Then, (gS]) follows from and (imj) . 

□ 

The scattering operator is defined as 

S :=]¥;]¥_. (4.13) 

COROLLARY 4.2. 

S = Po±. (4.14) 
Proof: This is immediate from fl4.3p because U* U = I. 

□ 

Let us denote by S± the restriction of S to 7^o±- 'S'_l is the physically relevant scattering 
operator that acts in the Hilbert space 7io± of finite energy fields that satisfy equations fl2.2p . 
We designate by I± the identity operator on Ho±. We have that, 

COROLLARY 4.3. 

S± = I±. (4.15) 

Proof: This follows from Corollary 14 .21 

□ 

The fact that Sj_ is the identity operator on 7io± means that there is cloaking for all 
frequencies. Suppose that for very negative times we are given an incoming wave packet 
g-*t^O(^_^ with G 7io±- Then, for large positive times the outgoing wave packet is given 
by e~**^°v9+ with y9+ = S±ip_. But, as S± = I, we have that ip+ = and then, 
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Since the incoming and the outgoing wave packets are the same there is no way to detect 
the cloaked objects K from scattering experiments performed in VL. 



5 The Proofs of Theorems 2.3 and 2.5 



It follows from (12.81) that the transformation matrix (12.121) is given by, 

= + - '-^) (y - c.)^(y - c,)^ 

yeBb^{cj),l<j<N, (5-1) 



_ ^ 113)3 \ I |Af 

The determinant is equal to 



A(y) = 9',i\y - c,|) y e B.^ic,), l<j<N. 



(5.2) 



A{y) = l,yeRl\uf^,B,^{c,). 

This result is easily obtained rotating into a coordinate system such that, y — Cj = (|y — 
Cj|,0,0) [n]. Then, by (l2:T8l) the matrices e^'',fi^'' are degenerate at dK and by (12191) the 
matrices ex,^, fixu are singular at dK. Moreover, by (12. 16^ 15. 15.21) there is a constant C such 
that. 



|e (x)| < C, |/''(x)| <C,^en. (5.3) 

Proof of theorem \2.3[ Let us denote, 

Q := {(E„Ho) e Ho : (Eo,Ho) = f/*(E,H)for some (E, H) G C',{Q)} , 

where U is defined in ( HAlh . Let us prove that Q C D{Ao). Note that Q C C^Qq \ 
Uf=i95b^(cj)). Suppose that (Eo,Ho) G Q. Then, 

«o( Ho ) eC(^]o\uf=l95,^(c,)), 
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where the derivatives are defined in classical sense on Qq \ UjLidBi,.{cj). Furthermore, by 
the invariance of Maxwell's equations, 

^° ( Ho ) = ^ ( Ho ) y ^ ^° \ ^Um, (c,). (5.4) 
This implies that. 



ao ( S° ) e no. (5.5) 



Ho 

But, as (Eo, Ho) have continuous tangential components at UjLidBh^icj) the function in the 
left-hand side of (15. 4p defined for y G Qo\^jLidBb^{cj) with the derivatives in classical sense 
actually coincides with the distribution that is obtained when the derivatives are taken in 
distribution sense in Rq. Then, by (15.51) we have that (Eo,Ho) G D{Ao), and 



/I / Eo \ / Eo 

^0 XT = cto 



HI ""U 1 TT 

/ V ^0 

with the right-hand side defined as indicated above. By (15. 4p this implies that, 

Ao\Q = U*anU, (5.6) 

or what is equivalent, that 

an = UAo\QU*. (5.7) 

It follows that to prove that oq is essentially self-adjoint and that (I2.42p holds we have to 
prove that Aq is essentially self-adjoint in Q. 

In the proof below we denote by an the formal differential operator. 

Suppose that (Eq, Ho) G Cq{Qo). Then, 

h:)-''-(h)-*"(h)-<h:)- (5*^) 

Hence, arguing as above, but in the opposite direction, we prove that, 

E \ J J ( Eo 



an\ H J = ^""M Ho ' ^ ^^'^^ 
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where the function in the left-hand side with the derivatives taken in classical sense in 
fl \ U^^^dBb_.{cj) actually coincides with the distribution obtained when the derivatives are 
taken in distribution sense in Q. Let us now introduce a mollifier. Take / G C^(]R^), /(x) = 
1, |x| < 1/2, /(x) = 0, |x| > 1, / f{x)dx = 1, and define /„(x) = re" f{nx). Denote, 

(h) :=^^(h)(^-^)^"(^)^^^^o°°^^)' ^5-10) 
for n large enough, where we have extended (E, H)^ to by zero. Moreover, 

(H)=-i™(H)„- <^-"' 

and 



s- lim an ( ^ j =s-\im^ f ^ j (x - z)/„(z) rfz = f ^ ) (x). (5.12) 



n^oo 



The limits in (15.111 15.121) are in the strong convergence in Tin. It follows from (15.61 15.81 
15.111 15.121) that (Eq, Hq)"^ can be approximated in the graph norm of Aq by functions in 
Q. Then, it is enough to prove that Aq is essentially self-adjoint in Cq{Qo). But as Aq is 
essentially self-adjoint in C^(Mo) it is also essentially self-adjoint in Cq(Mq). Hence, it is 
enough to prove that any function in Cq(Mq) can be approximated in the graph norm of 
Aq by functions in Cq{Qq). To prove this take any continuously differentiable real- valued 
function, /, defined on Mq such that, f{y) = 0, \y\ < 1 and f{y) = 1, \y\ > 2. Then, for any 



we have that. 



and moreover, 



flf{n\y-c,\) (^l)eCl{Qo), 

7 = 1 V / 



.,i.n/H^-.i)(H:)-(H; 
''»ri-f("iy-"='i) ( H„ ) ="» ( H 

.7 = 1 ^ ^ ^ 



5 
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where by s-lim we designate the strong hmit in Tig- This completes the proof that with 
domain Cq{Q) is essentially self-adjoint and that (12.421) holds. 

The unitary equivalence given by fl2.42p implies that Aq has the same spectral prop- 
erties that Aq. Namely, it has no singular-continuous spectrum, the absolutely-continuous 
spectrum is R and the only eigenvalue is zero and it has infinite multiplicity. Moreover, 



nn± := (kerneMo)" = | h J ^ ' d^x " ^ °' d^x'^ " ^ J ' ^^'^^^ 

□ 

Proof of Theorem \2.5[ - Let us denote by a the closure of a, with similar notation for a^, aj,j = 
I,-- - ,N. Then, 

a = Aq a], 

where we used the fact that as is essentially self-adjoint, = Aq. The adjoint of a is 
given by. 



and 
for 

( H. ) ( hO ^ ^'-''^ 

where the derivatives are taken in distribution sense. 

Let us denote by ]Cn± := kernel(2 =F '■— kernel(i =F «j) the deficiency subspaces of 

qq and aj,j = 1, ■ ■ ■ , A^. Since is essentially self-adjoint /Cn± = {0}. Let )C± := ©jLi^ji 
be the deficiency subspaces oi '■= ®^=iCLj- Suppose that /C± have the same dimension. 
Then, it follows from Corollary 1 in page 141 of [19j that there is a one-to-one correspondence 
between self-adjoint extensions of and unitary maps from /C+ into /C_. If V is such a 
unitary map, then, the corresponding self-adjoint extension Akv is given by, 

D{Akv) = {(p + (p+ + V(p+ : if e D(a^), ip+ e /C+} , 
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and 

Akv^ = ok^ + - iV(p+. 

Hence, since }Cn± = {0} and a = A^(BaK there is a one-to-one correspondence between self- 
adjoint extensions of a and unitary maps, V, from /C+ into /C_. The self-adjoint extension 
Ay corresponding to V is given by. 

Ay = An® Akv- 

6 Generalizations 

In this paper we considered high order transformation media obtained from singular transfor- 
mations that blow up a finite number of points, by simplicity, and since this is the situation 
in the applications. Suppose that we have a transformation that is singular in a set of points 
that we call M and denote now Vt^ := Mq \ M. What we really used in the proofs is that 
W^'2(M3) = Wo'^(no) where Wo'^(l]o) denotes the completion of C^iVLo) in the norm of 
W^'^(Mq). We also assumed that Eq'^^^q'^ are constant. What was actually needed is that 
Oq is essentially self-adjoint. Our results hold under this more general conditions provided 
that in fl4.2[ 14. 3p and (14.141) we replace Po± by the projector onto the absolutely-continuous 
subspace of Aq and that we assume that D{Aq) fl T^oac C W^'^(]Rq), where we have denoted 
the absolutely-continuous subspace of Aq by TYoac- Moreover, S±_ has to be defined as the 
restriction of 5* to Tioac and in fl4.15p has to be the identity operator on Tioac- Note that 
under these general assumptions Aq could have non-zero eigenvalues and singular-continuous 
spectrum. 

For example, W1'2(M3) = W]^'^{VLq) if M has zero Sobolev one capacity ^23 EHl ED- 
Moreover, assume that the permittivity and the permeability tensor densities e^^^^q^ are 
bounded below and above. Under this condition is essentially self-adjoint. Furthermore, 
let us denote by Tt^q the Hilbert space of finite energy solutions defined as in (12.31 1) but 
with = fiQ^ = S^'^. Let v4o,7tIo± be, respectively, the electromagnetic generator in Hq 
and the orthogonal complement of its kernel. We have that Tio and TtIq are the same set 
of functions with equivalent norms. Furthermore, D{Aq) = D [Aq), kernel Aq = kernel Aq. 
Moreover, (Eq, Hq)^ G 7^o± if and only if Eqa = EJ=o ^oa^Eom^ Hqa = Y.l=o /^oamHom for 
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some (Eo,Ho) G Hq^. As [211 [22] ^(^o) n C W1'2(R3) we have that D{Ao) n Hq^ C 
W^'^(]Rq) if Eq, fiQ are bounded operators on W^'^(]Rq) and this is true if the derivatives 
■jf-Eo, jf-fiQ are bounded operators on Hq for p = 1,2,3. Note, furthermore, that Hoac C 

Hq±. 

7 Conclusion and Outlook 

We gave a rigorous mathematical proof, in the time and frequency domains, that first and 
high order electromagnetic invisibility cloaks actually cloak passive and active devices in a 
very strong sense. This puts the theory of cloaking in exact transformation media in a firm 
mathematical basis that will allow us, in the next step forward, to analyze the stability of 
cloaking in the approximate transformation media that are used in the applications. 

A novel aspect of our results is that, as we prove cloaking for transformations media 
that are obtained from general anisotropic materials, we prove that it is possible to cloak 
objects contained in general anisotropic media. A very important case of anisotropic media 
are crystals. Suppose that we wish to cloak an object that is contained inside a general 
anisotropic medium, or a crystal, with permittivity and permeability tensors, respectively, 
^o'^ and /ig^. To cloak the object we proceed as follows. We have to coat it with a meta- 
material whose permittivity and permeability tensors are obtained using the transformation 
formula fl2.16p -for both of them- putting in the right-hand side, respectively, the permittiv- 
ity and the permeability of the general anisotropic medium, or the crystal. Our theory shows 
that this object with passive and active devices will be cloaked inside the general anisotropic 
medium, or the crystal. As we already mentioned, this is a new result, that shows that it 
is not necessary to transform from isotropic media. It is possible to transform from general 
anisotropic media. This opens the way to many interesting potential applications, not only 
for cloaking, but also for guiding electromagnetic waves under quite general circumstances. 
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